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ABSTRACT 


When low frequency underwater sound interacts with the 
Arctic ice cover, not only will it be reflected from the plane and 
scattered in all directions from roughness elements, but it will 
also be diffracted at leads and reradiated from flexural waves in 
the ice. These phenomena have been studied in an anechoic tank 
by pulse transmission from an underwater point source to a 
series of large floating acrylic plate models, each representing a 
different type of ice cover. The flexural wave speed, the plate 
and lead dimensions and the acoustic roughness are accurately 
scaled, and the specific acoustic impedance contrast 1s 
approximately modeled by the selection of the acrylic material. 
The physical contributors to the gross reflection coefficient and 
backscattering strength are identified and compared for models of 
a plane ice layer, an Arctic ice pressure ridge, edges of leads, 


and a rubble field of ice. 
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[. INTRODUCTION 


A. PROPAGATION OF SOUND IN THE ARCTIC 

In Arctic waters the existence of a positive sound speed 
gradient in the water column which includes a strongly positive 
near-surface portion which causes a shallow sound duct to be 
formed, forces sound energy to be repeatedly refracted upwards 
to the ice surface. Figure 1.1 illustrates a typical ray diagram 
for sound transmission in the Arctic. Long-range propagation 
measurements in the Arctic Ocean show that the low-frequency 
attenuation and backscattering strengths are more than an 
order of magnitude greater than in the open ocean (Mellen, 
1986). In order to predict the range at which a nolse source 
may be detected in the ice-covered Arctic, it is first necessary 
to understand what happens to the incident sound when it | 
interacts with the highly variable ice cover. To do this requires 
a basic understanding of the acoustical characteristics of the ice 


cover. 


B. DESCRIPTION OF TRE VARCTIGwiCh ee is 


In general, the Arctic ice surface is variable in thickness, 
density, elastic constants and shape. As illustrated in Figure 
1.2, the sea ice can be partitioned into the following acoustically 


significant categories: level ice (70 %), pressure ridges and 
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IMTERACTION OF LOW FREQUENCY SOUND WITH THE ARCTIC ICE CANOPY 
POTENTIAL CONTRISUTORS 
FLEXURAL WAVES OW ICE 


AIR ICE SAIL 


— 
An 7 -™ LEAD 





WW B—> “ply os 


RUBBLE FIELD 


WATER oN wo 


K K Ss Ss 


O = Diffraction from leading and trailing edges of plates 
K = Koel Scatter from deterministic ridges (ice keels) 


S = Scatter from steep-sioped. randomly-rougn, rubdle frelds 


6 = Boundary wave at steeo-sloped, randomiy-rougn, rubble fields, near grazing forward scatter 
Subsonic 


F =Flexural (surface) vrave, subsonic, evanescent in water, radiating into the sirsosce 
above the plate 


Figure 1.2. Scattering of low frequency sound from the 
Arctic ice canopy. Acoustically significant features of 
the ice canopy include level ice, pressure ridge keels, 
hummock fields and edges of leads and polynyas. 
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rubble fields (25 %), and leads (5 %). The number in 
parentheses is the approximate percentage of the Arctic ice 
cover each category represents (Weeks et al, 1971). This 
experiment is concerned only with kh < 1, where h is the ice 
thickness and k is the hcerizontal wavenumber. Each of these 
ice categories is briefly described below. 

1. Level Ice 

Sea ice 1s termed level if its surface changes less than 
1 part in 40 (Wadhams et al, 1985, p.1072). The thickness of 
level ice depends on the age of the ice, with first year ice 
generally less that 2 m thick and multi-year ice having a mean 
ice thickness of 2.5 to 3 m (Ackley et al., 1976, p. 9). 

2. Pressure Ridges 

Ridges are formed by the interaction of adjacent ice 
floes. The type and size of a ridge is a function of the thickness 
of the interacting ice and the degree of compressive or shear 
forces being exerted on the ice floes, 

Hibler et al (1972) compiled above-ice data from visual 
and laser profilometry and underice data from submarine 
underice profiles to develop statistical models of pressure ridge 
keel depths and the spacings between ridges. To insure that 
only pressure ridges are dealt with in the above ice and 
underice profiles, a cutoff height (or depth) was used, i.e., 
ridge height and keel depth distribution data are reported from 
a set that includes all ridges greater than a given cutoff. The 


average pressure ridge keel depth reported by Hibler et al 


iS 


(1972) from the set of all ridges greater than 6.1 m, was 

9.6 m in the Central Arctic and 11.4 m in the Canadian 
Archipeligo. They further report that the ratio of sail height to 
keel depth is on the order of 1:5. They report the average 
spatial distribution of pressure ridges as 4.3 ridges/km in the 
Central Arctic and 9.6 ridges/km in the Canadian Archipeligo. 

Hibler et al (1974) reported that ridging intensity is 
seasonally dependent. They determined that the spatial 
distribution of ridges in the Beaufort and Chukchi Sea regions 
was 1.83 + 0.6 ridges/km in March, and 0.99 + 0.61 ridges/km 
in October and November. Furthermore, they observed that 
all the ridges investigated had simular slope angles, such that 
the ratio of the base width to the ridge height was 
approximately four. 

A specific multi-year pressure ridge was the object of 
investigation by Kovacs et al (1973). This particular ridge had 
a sail height of 4 m and a keel depth of 13 m. The multi-year 
pressure ridge differs from the first year pressure ridge in that 
the latter is poorly bonded and contains many spaces between 
the blocks of ice that comprise the ridge. In a multi-year ridge 
these voids eventually are filled with ice due to the seasonal 
melting and re-freezing of the ice. 

Wright et al (1978) report values for the keel depth 
to sail height ratio and keel width to keel depth ratios slightly 
different from those above, i.e., 3.2:1 for keel to sail ratio and 


3.3:1 for the width to keel depth ratio. They also report that 


16 


keel depth statistics show little or no seasonal variability 
probably due to the fact that the underice surface does not feel 
the effects of the ablation of the ice surface during the warming 
months. 

In the present model the ice ridge is assumed to be a 
triangular cross sectional (prismatic) wedge below the plate and 
a similar smaller prism above the plate. The fine structure is 
omitted in this low frequency model. 

3. Hummock Fields 

In areas of severe shear forces ridges can be pushed 
together to form a series of adjacent ridges which would appear 
from above like a farmer's plowed field (Hibler and Ackley, 
1975). This type of non-random rubble field occurs where the 
forces acting on the ice are directional such as in a narrow 
strait or butted against a land mass. In another situation 
where the forces are less constrictive, such as in the Central 
Arctic Basin the ridges could be pushed together to form a 
nearly random field of jumbled blocks interspersed with large, 
nearly vertical ridge segments. It is this latter description 
which was modelled for this experiment. 

4. Leads and Polynas 

An open lead or polynya in the ice cover is an area 
where there is no ice. A covered lead exists where the 
thickness of the ice 1s much less than the that of the 
surrounding ice floes, 1.e., below 1 m thickness 


(Wadhams and Horne,1980). Leads are generally considered to 


shy 


be cracks in the ice where the length of the crack 1s much 
greater than the width. Polynyas on the other hand are large 
open areas up to several kilometers wide. For this experiment, 
open leads of 6, 12 and 18 m width have been modeled. A 
single edge of an open polynya has been modeled by using a 
single plate edge. Chapter III describes these models in greater 


detail. 
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[ee Reon, Cee eerie AND EQUIPMENT 


A. OCEAN PHYSICS LABORATORY AND ANECHOIC TANK 

The anechoic tank provided a low noise environment 
sufficient to conduct this experiment. The dimensions of the 
anechoic tank are: length - 7.3m, width - 1.8 m, depth - 
2.2m. Aluminum-impregnated rubber conical sound absorbers 
line each wall. Pulsed signals were transmitted and the signal 
return was identified, isolated and analyzed utilizing a judicious 
choice of a sampling window at the signal processing end. 
Figure 2.4 illustrates the equipment set-up utilized throughout 
the experiment. An equipment list Is provided in Table 1. All 
equipment used is available on the commercial market with one 
exception — the Random Surface Measuring Device (RSM). This 
device, which was designed by a former NPS student LT. S. 
Kasputis (1984), allows the accumulation of a large data base 
for statistics on a given surface. It can determine heights on 
slopes of less than 70 degrees with a precision of 20 um. 
Additional equipments which merit a more detailed description 


are discussed below. 


ey 


HP-467A 
Pwr Amp 


ITHACO 1201 
Low Noise Pre-Amp 


. HP-3314A 2 0 C 
[| ee Function Generator 
D > Oo 


0000000000 00000 


Interface Technology 
Timing Simulator /Word Generator 


IBM PC/XT 


A/D Converter 


sampling frequency 
= 1 MHz 





Figure 2.1: Schematic diagram of the signal processing 
flow path. 


20 


B. TRANSMIT AND RECEIVE SYSTEMS 
1. Hewlett Packard Function Generator, Model 3314A 





The 3314A Function Generator is a digital, multi-mode, 
programmable function generator which features preset Sine, 
Triangle and Squarewave functions and programmable arbitrary 
waveforms from 1mHz to 19.999MHz and a maximum output 
voltage of 10 volts peak-to-peak. The function generator was 
utilized in the N CYCLE MODE which outputs a burst of N 
complete cycles, starting when an internal or external trigger. 
edge is received. The excellent stability of the function - 
generator permitted repeated data runs in minimal time. 

2. Nicolet FFT Wave Analyzer , Model 660B 

The Nicolet wave analyzer is programmable and performs 
its own analog to digital conversion. The sampling rate is 
limited to 250,000 samples per. second. It has a built-in CRT 
screen that allows real time observation of either the received 
time signal or its frequency spectrum. This equipment was 
utilized in parallel to the signal processing equipment to permit 
initial detection of the signal of interest. 

3. Interface Technology Timing Simulator, Model RS5-648 

The timing simulator was used to trigger the function 
generator, the Nicolet and the signal processing equipment. It 


has a 50 nanosecond timing resolution. 


zu 


4. Celesco LC-10 Hydrophone 

LC-10 hydrophones were used as both transmitter and 
receiver. The specifications for the experiment required a point 
source (ka « 1, where ‘k' is the wave number and ‘a’ is the 
linear dimension of the transducer) and receiver with rise and 
decay times sufficiently short to permit the reception of 
undistorted pulsed signals in an environment where the time 
window was critical. The time duration of the pulse is 
important since all of the desired signal must be sampled in a 
time window free of the unwanted arrivals from the tank or 
plate reflections and diffractions from the edges of the model. 
Additional requirements were a relatively high and flat 
free-field voltage sensitivity and an omnidirectional directivity 
pattern in the frequency range of 30 kHz to 100 kHz. The 
LC-10's directivity pattern was omni-directional to within + 2 
dB. The greatest disadvantage of the LC-10 as a source is its 
low output power, however, the combination of short ranges 
(less that 2 m) and the low ambient noise level in the anechoic 


tank permitted its use with excellent results. 


C. SIGNAL PROCESSING 2OUIE Mii 


Signal processing was accomplished using an IBM PC/XT 
personal computer fitted with a Computerscope Model ISC-16 
signal Processing system. The Computerscope [SC-16 system is 
a fully integrated hardware and software package designed to 


permit the IBM to perform as a data acquisition and analysis 


oe 


laboratory instrument. The [SC-16 package consists of a 
16-channel analog to digital converter board, an external 
Instrument Interface and Scope Driver software. The package 1s 
capable of receiving 16 channels of data input at an aggregate 
sampling rate of up to 1MHz. Digital conversion is achieved 
With 12 bit accuracy over an input range of -10 to +10 volts. 
The Scope Driver software permits the IBM CRT to emulate a 
digital storage oscilloscope with numerous waveform 
manipulation and display modes, and to output high resolution 
graphics to a dot matrix printer or X/Y plotter. Figure 2.2 is a 
graphic produced by the Computerscope system and illustrates 
the systems excellent signal processing capabilities. The package 


1$ capable of directly accessing 640K of RAM. 


US 


TABLE | 


EQUIPMENT LIST 


Abbreviation 


Timing Simulator 


Function Generator 
Power Amplifier 


Source or Receiver 


Pre-Amplifier 
Nicolet 
Voltage Amplifier 


Signal Processor 


Scope 


Nomenclature 


Interface Technology Timing 
sniula ory a Generator, Model 


Hewlett Packard Function Generator, 
Model HP-3314A 


Hewlett Packard Power Amplifier, 
Model HP-467A 


Celesco LC-10 Hydrophone 


Ithaco Low Noise Pre-Ampifier, 
Model 1201 


Nicolet Dual Channel FFT Analyzer, 
Model 660B 


Hewlett Packard Voltage Amplifier, 
Model HP-465A 


IBM PC/XT with Computerscope 
ISC-16 Data Acquisition and Analysis 
Package 


Kikusul Oscilloscope, Model COS 5060 
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Time (msec) 
Figure 2.2: 
[Mimjiewcaonmaain Udltage response to an input sinusoidal 


ues sarivemeat /.U0/ volts RIMS.) Sintisoid frequency 
reoo ke 


ioe ose. u(t) = A sintZart) [ilt) — 1t-277)] 


Peak Identification: 
& = bireer signal 
— Surface reflected signal 
(3) — Specular scatter from sides of tank 
(4) - Diffraction from corner of water-air interface 
and side of tank 
© — Specular scatter from the bottom of tank. 
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Il. DESCRIPHON OF PHY SIGaL Onis 


Before construction of the “ice” models could begin, it was 
necessary to select the material which best approximated the 
physical characteristics of Arctic ice. The internal structure of 
sea ice, such as its saline content and the distribution of melt 
inclusions, is beyond the scope of our interest; furthermore for 
the low frequency regime which is our interest, average values 
provide a satisfactory acoustical description. Therefore, for this 
experiment it was necessary to emulate the bulk properties of 
the ice cover. Acrylic plate was chosen because it has several 
important physical properties that parallel closely those of Arctic 
sea ice. Specifically, the longitudinal wave speed and flexural 
wave speed (for scaled thicknesses), Poisson Ratio, and the 
characteristic impedance (pc) for acrylic plates all fall well 
within the measured range for Arctic sea ice. The density of 
the acrylic plate is about 25% too high and the shear and 
compressional attenuations of the plates are somewhat higher as 
well. Table Il summarizes the comparison between Arctic ice 
conditions and our laboratory acrylic “ice sheet” floating in a 


fresh water tank. 
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TABLE I] 


COMPARISON OF THE PHYSICAL PROPERTIES OF 
ARCTIC ICE TO THE ACRYLIC MODEL LABORATORY 


PROPERTIES 


CHARACTERIST! 


fer OUNGS MOD 
E, (Pa x 107) 


Zee DENSITY 
p> (kg/m?) 


5. POISSON'S RATIO 


4. BULK COMPRES- 
SIONAL WAVE 
SPEED (m/s) 


5. COMP. WAVE 
ATTENUATION (a) 
(dB /m), F(kKHz) 


6. BULK SHEAR 
WAVE SPEED (m/s) 


7. SHEAR WAVE 
ATTENUATION (c) 
(dB/m), F(KHz) 


CONDITIONS 


ARCTIC 


RANGE TYPICAL 


ee ee ee eo eee oe = ee ee 


7 7-0 4.1 


(1) (2) 
760-960 930 
(4) (4) 
0.30-0.36 0.33 
(5) (5) 
2800-3500 3150 
(6) (6) 
0.3F-0.05F 0.06F 
(20) (20) 
1500-1860 1700 
(6) (6) 
1.8F-0.30F 0.36F 
(20) (20) 


Jat 


LAB 


1148 
(3) 


0.316 
(3) 


2240 
(3) 


0.424F 
(3) 


1373 
(3) 


0.853F 
(3) 


MODEL/ICE 
RATIO 


0.46-2.41 


Uo | 


OF oo 1eO> 


URDU > 1 ea5 


Y 41sSes. 45 


0.74 - 0.92 


Oa ee 2) 


1) 


ie 


Ze 


13. 


TABLE II (CONT'D) 


COMPARISON OF THE PHYSICAL PROPERTIES OF 
ARCTIC ICE TO THE ACRYLIC MODEL LABORATORY 


BULK CHAR. IMP. 
202, 
(kg/m?-s) X 108 


FLEXURAL WAVE 
SPEED (m/s) 


FLEXURAL WAVE 
ATTENUATION 
(dB/m), F(KHz) 


PLATE THICKNESS 
h (m) 


THICK/WAVELENGTH 0.007-0. 1 


h/A 


VEL. RATIO 
Co/Cy 
(ICE/WATER) 


CONDITIONS 
1.06-3.54 2.07 
(8) (9) 
252 =255/ 898-1270 
(11) (12) 
Data not 
available 
0.5-6.0 276 
(4,7) (7) 

0.09 
(14) (15) 
0.97-2.48 1.54 
(17) (18) 
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ane) 
(10) 


700-1150 . 
(3) 


1.5F 
(3) 


0.00318 
(3) 


0.04-0.11 
(16) 


1.54 
(19) 


O./357532 


0.49-2.78 


1: 600 


O4= ee 


Oro l—ioe 


TABLE I] (CONT'D) 
COME ole Ot semen cCnl PROPERTIES OF 
PRCM CCE TO Tir ACR YEIC M@BEL LABORATORY 


CONDITIONS 
WATER CHARACTERISTICS 
14. TEMPERATURE (-1.6) to (-1.8) 20 
ce (1) (3) 
15. SOUND SPEED 1440-1485 1440 14814 1.00-1,03 
C,,(m/s) (22) an (3) 
16. DENSITY 1000-1040 1026 1000 0.96-1.00 
p,, (kg/m?) (21) (21) (3) 
FREQUENCY RANGE 
OF INTEREST 0.05-0.14 30-80 600: 14 


F, (KHz) 


LIST OF REFERENCES FOR TABLE II (the complete reference for all 
of the below references is in the LIST OF REFERENCES). 


Weeks, 1976. 

Typical value. 

Measured laboratory value. . 

Pemeyeet al, CRREL Report 7o—-id, 1976. 
ochwartz and Weeks, 1977. 

Hunkins, 1960. 

Wadhams et al, 1985. 
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TABLE II (CONT'D) 


COMPARISON OF THE PHY signer OPE Ritoned 
ARCTIC ICE TO THE ACRYLIC MODEL LABORATORY 


Calculated 


Calculated 
Calculated 
Calculated 


Calculated 
Calculated 
Calculated 
Calculated 
Calculated 
Calculated 


Calculated 
Calculated 


CONDITIONS 


using: p9 = 760 to 960 kg/m’, 

and cp = 1395 to 3688 m/s. 
using: p9 = 930 kg/m%and cy = 2224 m/s. 
using: p9 = 1148 kg/m%and cy = 2240 m/s. 


using: v = (2nfh) 2B (m/s), where: 
hmin = 0.5 m, hmay = 6.0 m 


fmin = 20 Hz, fmax = 140 Hz 
vee Cues Bmax = 392.44 where 
B is the flexural wave constant (m/s) ”: 
using h = 2.0 m and B = 25.34 (m/s) ?: 
using h = 3.18e-3 m and B = 23.98 (m/s) 2 
using h = 2.0 m and f = 50 to 140 Hz. 
using h = 2.0 m and f = 95 Hz. 
using h = 3.18e-3 m and f = 30 to 80 kHz. 
USING Cyan meio o, Cmax = 3688 (m/s) 
Cmin = 1440, Cee cer noe 
Using seu=— dao. =) 2224S) 
using c = 1481, c = 2240 (m/s). 


McCammon and McDaniel, 1985. 
Myers et al, 1969. 
NAVOCEANO TECHNOTE, 1974. 


A. MODELS OF SMOOTH ICE AND ARCTIC OPEN LEADS 

scaling of the parameters to be measured such as 
diffraction strength, coherent forward scattering and flexural 
wave propagation for a particular frequency and ice thickness 
range in the Arctic can be adequately achieved. This requires 
that the physical properties of the model be similar to those 
Same properties in the Arctic. Once this is accomplished, what 
is required is that the scaling of the ratio of model length to 
laboratory wavelength be the same as Arctic conditions. As 
previously stated, the acoustically significant physical properties 
of Arctic ice are well matched in acrylic. Therefore, modeling 
of the arctic for several of the phenomena studied here requires 
that the ratio of the acoustic wave length (A), and the ice 
thickness (h), be the same for the two situations. Medwin et 
al. (1984) successfully demonstrated the accuracy of this . 
modelling technique. The frequency range of interest for sound 


propagation in the Arctic for this experiment is 50 to 150 Hz. 
\ = speed of sound / frequency 
Lh/Al ice = Lh/A] acrylic 


i acrylic ~ Lh] acrylic / Lh/Alice 


al 


Assuming a typical Arctic level ice thickness of about 2 m, then 
1/8 in acrylic (3.125 mm) corresponds to a model scale of 


about 600:1. 


Therefore, 
when f jce = 50 Hz, f acrylic = 50 kHz, and 
when f jeg = 150 Hz, f acrylic = 90 kHz. 


Since the density of acrylic is slightly greater than water 
it was necessary to devise some means of “floating” the plate in 
the tank without suppressing the generation of flexural modes in 
the plate. The solution to this problem was actually quite 
simple. A 1.27 cm high bead of silicon sealer laid around each 
edge of the plate effectively turned the plate into a 
flat-bottomed boat which floated freely on the water. The 
dimensions of the plate were 2.4 m by 1.5 m. These large 
dimensions helped to ensure that extraneous edge diffractions 
would not interfere with the signal of interest. To further 
minimize unwanted edge diffraction three of the four edges were 
tapered. On the fourth edge (along the width) was glued a 10 
mm high by 3 mm thick acrylic bar. It was this edge, when 
separated from a similar edge on a second plate by the desired 


distance, which was utilized to simulate an open lead. 
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PmeeviODEL OF A PRESSURE RIDGE 

A cross sectional view of the single pressure ridge model 
showing the significant dimensions and scaling is given in Figure 
3.1. The scaled dimensions of the sail and keel were designed 
to fit those of a typical multi-year pressure ridge as given in 
section I.B. The keel was manufactured by glueing eaeiner 
two flat strips of 3/8 in acrylic and machining to the specified 
dimensions and shape. The sail was machined from a single 
3/8 in thick acrylic strip. Both the sail and the keel were then 
glued using an acrylic solvent (K-Lux® solvent cement for 
acrylic) to a smooth acrylic sheet. Because the Arctic 
frequency range of interest (50 - 150 Hz) has wavelengths 
greater than 10 m, which is much greater than the Arctic sea 
ice thickness, it was not necessary to introduce micro-roughness 


elements on the laboratory model of the ridge. 


ee vODEL OF A RUBBLE FIELD 


The construction of the rubble field was the most difficult to 
achieve of the three models. Statistics of a rubble field, per se, 
were not available. However, a study of available literature on 
pressure ridges indicated that a typical mean ice plate thickness 
of 4.0 m with an underwater depth standard deviation of 
1.0 m are representative values for pressure ridges (Wadhams 
et al, 1985). A rubble field can be described as a field of 


pressure ridges which are randomly oriented with close (several 


SS 





Figure 3.1: Cross section of the pressure ridge model 
used to measure forward and back scatter from a ridge 
keel. Dimensions in parenthesis are scaled to the Arctic 
model. All materials are acrylic. 
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meters) but variable spacing (Bourke, 1986). Thus, the 
dimensions as stated above were used as guidelines for the 
construction of the rubble field. Using the scale 600:1 this 
corresponds to lab dimensions of 6.7 mm mean thickness and 
1.7 mm standard deviation. 

The rubble field was constructed from acrylic plates of 3 
mm and 9.5 mm thickness which were cut into strips ranging 
in width between 3 mm and 9 mm. The majority of these 
strips were then cut again along the width in 3 to 9 mm size 
strips to form randomly sized blocks of acrylic, and the 
remaining strips were broken into jagged pieces by chiselling 
with a sharp-edged hammer. Enough pieces of acrylic were 
made this way to cover a 3 m® surface. The dimensions of the 
rubble field were 0.6 m by 1.6 m, which when scaled to the 
Arctic would represent a 360 by 960 m surface area. The 
pleces were fastened to the top and bottom of the plate using 
cyanoacrylate adhesive. The rubble was applied so that it 
would match as closely as possible the dimensions of 6.7 mm 
mean thickness and the 1.7 mm ice bottom roughness standard 
deviation stated above. 

The statistics for the plate underside roughness were 
determined using the Random Surface Measuring device (RSM). 
Six 18 cm long traverses of the plate were completed and a 


total of 3200 data points were collected at an incremental 


Jo 


distance of 0.332 mm and with a height resolution of 20 um. 
A typical profile of the surface is shown in Figure 3.2. 

The surface data were processed to determine the mean 
height and standard deviation with respect to the mean height 
of the plate underside roughness. The probability density 
functions of the surface heights and slopes were obtained and | 
are illustrated in Figures 3.3 and 3.4, respectively. The best 
fit for both of these curves is indeterminate, but it is not 
gaussian. The acrylic smooth plate thickness is 3.1 mm. The 
mean thickness of the bottom roughness glued onto the plate is 
3.8 mm. The standard deviation of the bottom roughness is 
2./6 mm. A plate topside roughness of approximately one-half 
the thickness of the bottom roughness was also glued onto the 
plate for purposes of mass loading but its dimensions were not 
measured. These dimensions represent an Arctic rubble field 
mean ice thickness of approximately 5 m and an ice bottom 


roughness standard deviation of 1.65 m. 
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Figure 3.3. Height histogram ofthe bottom Sumiaceroi anne 


roughened acrylic plate that was manufactured to model 


an Arctic rubble field. The height dimensions are presented 


in actual model dimensions (mm) and in scaled (600:1) 


Arctic dimensions (m). 
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Figure 3.4. Slope histogram of the bottom surface of a 


roughened acrylic plated used to model an Arctic rubble field. 
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ocattering is defined as the re-radiation of incident acoustic 
energy by discontinuities in the physical properties of the medium 
(Urick, 1983). These discontinuities can be manifested in water 
as air bubbles, a rough air—water interface, or the rough 
undersurface of the Arctic ice cover. For the purpose of this 
work scattering from the Arctic ice cover will be composed of the 


following acoustically significant effects (see Figure 1.2): 


1. Pressure waves diffracting from edges and propagating 
in all directions. 


2. scattering in the specular direction, i.e., the angle of 
reflection is equal to the angle of incidence. 


3. Conversion of the incident acoustic energy to shear and 
flexural modes, and formation of compressional modes 
by refraction, in the plate. This energy can be 


re-radiated back into the water below, or the airspace 
above, the plate. 


The theoretical aspects of each of these effects will be discussed at 


length in this section. 


A. DIFFRACTION 

Diffraction can best be visualized using Huygens’ principle. 
When an acoustic wave encounters an obstruction such as the 
edge of a plate, every point on that edge will be the site of an 
expanding spherical Huygen's wavelet. The acoustic energy 
re-radiated in all directions from the edge in these spherical 


wavelets 1s called diffraction. 


Gy, 


Biot and Tolstoy (1957), in applying the use of normal modes 
as generalized coordinates in Hilbert Space to mechanical media, 
developed an exact, closed form, expression for the diffracted 
pressure wave from a rigid wedge or corner. Medwin (1981) 
modified the source function used by Biot and Tolstoy to a delta 
function source that is turned on at time t= 0, and simplified 
the equation for the diffracted pressure. Figure 4.1 shows a 
wedge to illustrate the typical cylindrical coordinates used in the 


diffracted pressure wave equation given below: 


(-Spc/476,,,) {8} expl-my/9,,] 


, td, 
pit} = poems (4-1) 
0 Pt 
Where, 
es sint{7/6, }{1+6+8,)] 
1-2 expl-7y/6, casli/6, l2620,)] + expl-eny/e,)  '4>2 


cot ile eee 


Y = arc cosh 
Zine 


ae 


where p(t) is the diffracted pressure at time t, p is the density 
of the medium, 5 is the source strength, and c is the sound 
speed in the medium. The source coordinates are (r,,6,,0) and 


the receiver coordinates are (r,6,z). 6, is the angle of the 


W 


wedge measured in the fluid region. The term (n+6+6,) is 
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Figure 4.1: (a) Wedge geometry for equations 4-1 thru 
4-3. (b) Unfolded geometry. 
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written for simplicity; the curly bracket in Equation 4-1 consists 
of the sum of the four terms obtained by using the four possible 
combinations of the angles, +8 and +0,. 

Physically, the sound energy emanating from the point source 
diverges spherically, intersects the ridge, and is diffracted to the 
receiver with the earliest arrival occurring at the least time, 


(t,), given by: 
ty = [Cr + r,)? + 22} 2/c. (4-4) 


oubsequent arrivals are due to diffraction of the spherical wave 
as 1t propagates over the wedge from elements of the wedge on 
either side of where the least time path intersects the wedge. 
The maximum pressure occurs at the least time and falls off 
exponentially with increasing time. 

The diffracted pressure field from the acrylic model, or the 
Arctic ice cover, can be determined by the linear superposition 
(Novarini and Medwin, 1985) of Equation 4-1 applied to a series 
of wedges that compose a specific feature on the model, or ice 
cover. For example, the model of the ridge keel is composed of 
three distinct wedges labeled 1, 2, and 3 in Figure 4.2a. Figure 
4,.2b shows the result of applying Equation 4-1 to each of the 
wedges, with r = ry, and 8 = 9, and summing. 

Consider the ridge keel as a linear system whose output is the 


diffracted pressure, then the system diagram 1s 
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60 120 180 240 300 
Time (jisec) 





Figure 4.2: (a) Diagram showing the three Biot~Tolstoy 
wedges that make up the model of the ridge keel. 

(b) Time domain impulse response calculated using Biot- 
Tolstoy theory and superposition. Relative times of 
arrival for a backscattering geometry. 
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s(t) > = n(t) 


where h(t) is the impulse response of the filter and H(f) is the 
Fourier transform of h(t) and is called the transfer function. 
The diffracted pressure response, p.(t), for the pulsed sinusoid, 


y(t) (see Figure 2~2), is the convolution of y(t) with h(t); 


pS (C7 eee 


Using the convolution theorem, this becomes, 
P,(f) = Y(f)H@) 


in the frequency domain. The pressure response in the 
frequency domain relative to direct free-field transmission at a 
reference range R, Is simply the filter transfer function multiplied 


by the reference range: 


Po(f)_ _ Y(f)H(f) 


Preilf) = Y(f/Ro Y(f)/Ro 


= H(f)R, 


In this treatment, the diffracted pressure from Equation 4-1 is 


given by: 


p(t) = (eS/4n)&(t)*h(t), 


ae 


in the time domain, such that, 


P(f) = (eS/4n) H(f) 


in the frequency domain. Now the relative frequency domain 
pressure response for any given source can be calculated from the 


Fourier transform of Equation 4-1 by: 


4nR 
Cae) = eS : P(f) 


Although there is no complete analytical transform for Equation 
4-1, the transformation can be performed using digital 


techniques. Equation 4-4 then becomes: 


Pa) = SP) (4-5) 


where At is the sampling interval and is equal to the reciprocal 
of the sampling frequency used in the digitization of Equation 4-1 
(Brigham, 1974). Medwin et al (1982), using a discrete Huygens 
interpretation, obtained a finite value for p(t=1,) and showed 


that this was related to p(t=t,+At) by 
p(t=t,)/p(t=1,+At) = 1.366. (4-6) 


Biot and Tolstoy develped their theory using infinite 


dimensions for the wedge, and rigid boundary conditions. The 


ad 


wedges that make up the ridge keel and leads in the laboratory 
models (and in the Arctic) have finite dimensions and real, 
impedance boundaries. Novarini and Medwin (1985) report that 
the Biot-Tolstoy theory with discrete Fourier transformation can 
be applied to finite wedges as long as the wedge separations are 
greater than one-quarter wavelength. Experimental verification 
can be found in Medwin et al (1982) and excellent agreement 
between theory and experiment for a rigid wedge is given in 
Medwin (1981). To date the effect of non-rigid boundary 
conditions on the Biot-Tolstoy diffraction theory has not been 
investigated, but it is expected that a non-rigid wedge will give 
less diffraction than a rigid wedge, due to the conversion and 
refraction of incident acoustic energy to acoustic energy in the 


medium of the wedge. 


B. ROUGH SURFACE FORWARD SCATTER FROM A RUBBLE FIELD 
The use of the Helmholtz theorem with the Kirchhoff plane 
wave approximation was first applied by Eckart (1953) to 
calculate the acoustic field scattered by a Gaussian rough plane 
surface. It has been shown that the coherent component of the 
scattered signal can be obtained by averaging the signal over the 
probability density function (PDF) of the rough surface (Clay and 
Medwin, 1977). For specular scatter, when Scattered = 
Gincident» in addition to the coherent component,<p>, due to the 
statistically weighted height-—dependent phase shift, there is an 
incoherent component due to the spatially dependent phase shifts 
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as determined by the spatial correlation of the heights. That is, 
<p¢> = <p>4 + Var(p) (4-7) 


where <p¢> is the total scattered intensity and <p>¢ is the 
coherent component. For a low roughness surface, kocos; < 1, 
the coherent component dominates and the coherent scattered 
intensity is approximately equal to the total scattered energy, 
i.e., <p> ~ <p4>,. Here k = 2n/A, o = root mean square (rms) 
height of the roughness elements, and 6; = angle of incidence 
measured from the normal. . 

The received signal from sound specularly scattered from a 
rough surface is the sum of the scattered signals from all of the 
scattering elements of various heights ensonified by the incident 
plane wave, and arriving at the receiver in the selected time 
window. The probability of occurrence of a specific height of a 
scattering element between J and ( + dZ is wa(%)d(%), where Z 
is the height above the mean level and w,(%) is the PDF of the 
rough surface. The effect of ¢ on the signal is to alter the phase 
by -2k{€cos0;. Equation 4-8 derives the mean value of the 
coherent reflection coefficient ,<R->, by multiplying the effect of 
exp(—2ik@cos8;) by w,(%) and integrating over all values of J. 


(Clay and Medwin, 1977) 


4/ 


oO 


ge A, | expl-2ik £c0s6,] Wal Iz (4-8) 


-O0 


where Rg is the plane wave smooth surface pressure reflection 
coefficient. 

Assume the rough surface has a Gaussian PDF, Equation 4-9. 
Then, substitution of Equation 4-9 into Equation 4-8 and 
evaluation of the infinite integral results In Equation 4-10 which 
describes the magnitude of the coherent component of the 


specular scatter due to the roughness elements 
wg) = 071 (20) *exp(-C2/[202}) (4-9) 
<R,>/R, = exp(-2k*o%cos 6) = exp(-g /2) (4-10) 


where g is defined as a roughness parameter dependent on 
frequency (k = 2nf/c), incident angle (6), and rms height of the 
roughness elements (a). 

This theory assumes several conditions which are not met in 
the Arctic or in this experiment, namely: incident plane waves, 
low roughness, 1.e., @ < 1, and no steep-sloped roughness 
elements. Spherical waves will diverge as 1/range while 
propagating; this factor will be corrected later (Equation 5-3). 
The fact that g > 1 in some of our experiments will result in 


<R,*> greater than <R,>¢; that is Equations 4-8 and 4-10 do 
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not measure the total scattered energy. Finally Equation 4-8 1s 
derived for gently sloped roughness elements, i.e., less than 10° 
slopes, whereas the Arctic, and our model, has roughness 
elements with slopes as great as 90°. Therefore, for a roughness 
parameter greater than 1 there will be additional energy at the 
receiver due to the increasing significance of the incoherent signal 
component. Hence, we would expect a received signal <p4> 


larger than theory predicts for <p>. 


C. REFLECTION FROM A SMOOTH PLATE 

The calculation of the plane-wave smooth surface reflection 
coefficient, R,, for a lossy, multilayered medium such as Arctic 
ice 18 a complicated problem due to its complex and anisotropic 
nature. The lossy, multilayered smooth ice plane-wave reflection 
coefficient was recently studied by McCammon and McDaniel 
(1985). They extended the Thompson—-Haskell matrix method to 
treat a layered viscoelastic medium, 1.e., ice bounded by two 
fluid half spaces, water and air. In their computer model they 
incorporated the density, the bulk compressional speed and 
attenuation, and the bulk shear speed and attenuation constants 
of ice. They reported that ice thickness and shear attenuation 
were the critically important paramenters in determining Rg. 
They were kind enough to provide us with a copy of this program 
in order to derive theoretical predictions for the smooth acrylic 


plate reflection coefficient. 
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C. PLEAURAL WAVES Mi SIICO (rises 

Flexural waves in plates arise as a result of transverse 
vibrations at frequencies kh < 1, where k is the horizontal 
wavenumber and h is the plate thickness. There is a 
displacement perpendicular to the plane of the plate that 
propagates in a direction that is parallel to the plane of the plate. 

An excellent review of the derivation of the flexural wave 
equation for both uniform beams and plates is contained in Ross 
(1976). In his treatment he uses a non-dimensional shear 


parameter (T oS 
rT, = Y/[K,G(1-9°) ] (4-11) 


for the plate, which relates Young's Modulus (Y) to the shear 
modulus (G). The factor, Ka loa)e is included to account for 
Warping in the plate cross section resulting in a non-constant 


shear force. For a plate, 
Kp ~ 0.76(1 + 20/5). (4-12) 
where o is Poisson's ratio. 


Assuming that the spatial derivitive of thickness is 


negligible, the homogeneous bending equation for a plate is given 


by 


90 
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where y is the displacement perpendicular to the plane of the 


plate, c, 1s compressional wave speed in the plate, and h is plate 


p 
thickness. Two of the solutions to this equation are real, and 
two are imaginary. The general solution can be expressed as 


(Ross, 1976; Kinsler et al., 1982): 
y = [Ae Ikx + Belkx + Ce-iv% + Delyx]elot, (4-14) 


where w = 2nf, k is the flexural wavenumber = w/v;, vy; is the 
flexural wave phase speed, y = ik, where i = (-1)”, and 
A, B, C, and D are complex coefficients determined by the 
boundary conditions. 
1. Flexural Phase Speed 
Flexural waves travelling in a plate propagate at a phase 
speed given by w/k. Substituting the harmonic solution from 


Equation 4-14 back into Equation 4-13 gives: 


2 73 
Zi 24 Hae 2 ae 
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which .can be solved for the flexural wave phase speed, vy, by 
substituting w/ Vs for k, and solving the resulting quadratic 


equation. This gives 





f 2A 
= ae 
| = 2 
OO Cy l2 
(4-16) 
2 
- ne 
B= If, N55 
2 
é., ee 
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Note, that at low frequencies, A is dominated by 1/w*, therefore 
4AC dominates B* and (4AC)” dominates B resulting in the 


familiar equation for flexural wave speed given by Kinsler et al. 
(1982): 


GIC_N 
oe 


f vie (4-17) 





At low frequency the flexural wave phase speed is proportional to 
the square root of the frequency (see Figure 4.3). 
At very high frequencies, A > -(f,h?)/(c,°12), and the 


flexural wave phase speed becomes a constant, 1.e., asf > @, 
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Therefore the flexural wave phase speed is dispersive with (f)* at 
low frequencies, and non-dispersive at high frequencies. It is the 
proportionality for low frequency speed, v; ~ (wh)” ~ (h/d)®, 
that permits scale modeling to be feasible for flexural waves. 

2. Flexural Wave Group Speed 

Group speed v, is given by dw/dk. This can be obtained 
by solving Equation 4-15 for w*, then solving for the derivative 
of w with respect to the wavenumber k. This results in the 


following expression: 


Z 
2h 2 
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It can be shown that as w gets very small, 


Vg = dw/dk - 2[we h/¥12]" = 2vel, oman (4-20) 
As w gets very large, 
Vz = dw/dk > eaclils = Vel. large: (4-21) 
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That is, at low frequencies the group speed is twice the phase 
speed and at very high frequencies the group and phase speeds 
are the same. A plot of the flexural wave phase and group 
speeds as a function of frequency for a 3 mm thick acrylic plate 
in alr is given in Figure 4.3. 
3. The Effect of Entrained Mass on Flexural Wave Speeds 

The medium in which the plate is situated will affect the 
net radiation impedance seen by the vibrating plate. The 
reactance term of the radiation impedance can be viewed as the 
effect of the entrained mass of the medium directly surrounding 
the vibrating plate ( Kinsler et al., Chap 8, 1982). If the 
medium is air, then the effect of the entrained mass in the 
above equations will be negligable. However, if the medium is 
water, the entrained mass will be significant and a correction to 
the equations for phase and group speeds will have to be used. 

Letting e be the ratio of entrained mass to the plate 


mass, the terms in Equation 4-13 become modified as follows: 
oe ee and 
h? > h@/(1-e). 

The flexural wave phase speed can then be determined in the 


same manner as described above, with the resulting low 


frequency approximation given by: 
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vy? = weph/[12(1+e)]/ , kh < 4 C=2>) 


At high frequencies, the radiation reactance goes to Zero, 
therefore the effect of the entrained mass disappears and the high 
frequency approximation is the same as Equation 4-18. 

Entrained mass will reduce the flexural wave phase and 
group speed at low frequencies, but as frequency increases the 
effect of the entrained mass becomes less and becomes Zero as 
frequency approaches infinity. 

4. Boundary Laver Formation due to Flexural Waves. 

The condition necessary to achieve efficient mode 
conversion between compressional waves in the water and 
flexural waves in the plate is to match the horizontal component 
of the compressional wavenumber (x) in the water with the 
flexural wavenumber (k,) in the plate. The condition of 
matching horizontal wavenumbers is known as coincidence. Both 
in the laboratory model and in the Arctic, the speed of flexural 
waves is less than the compressional wave speed in the water 
underlying the plate for all frequencies. Therefore, the condition 
for efficient mode conversion is never satisfied due to the subsonic 
nature of the flexural wave compared to waterborne waves, 1.¢., 
ke is always greater than x. 

Any modes re-radiated into the water from a subsonic 
flexural wave in a plate will be evanescent, that is, the energy 
will decay exponentially with distance from the plate instead of 


propagating as a harmonic wave. In a recent paper Tolstoy 
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(1986) showed that this should form an acoustic boundary layer 
below the plate and developed expressions for the pressure field in 
this boundary layer for smooth and rough plates. For a point 
source located in the water at a depth z, below a smooth plate 
and a receiver at a depth z, the magnitude of the pressure in 
the boundary layer (neglecting attenuation) relative to freespace 
acoustic pressure 1s given by: 

ee expla (cten)}_ (4-23) 

ED ae kf V¢ vel ky ° 
where r is the horizontal distance between source and receiver, 
P, is twice the freespace acoustic pressure at r, k and c are the 
compressional wavenumber and wave speed In water, 
respectively, k;y is the flexural wavenumber for the plate, and 
g=(k,* - k?)”. If attenuation in the plate is significant, this 
equation should be multiplied by exp{-ar) where a is the 


flexural wave attenuation coefficient in Np/m for the plate. 
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V. “SsPECULAR SCA: 


A. EXPERIMENTAL PROCEDURE 
1. Data Collection 

Laboratory measurements of the smooth plane wave 
pressure reflection coefficient (Ry), the rough surface coherent 
pressure reflection coefficient (Ry), and the wedge forward 
specularly scattered pressure were made at specular angles in 
the water-filled tank described in Chapter I]. Both the source 
and receiver were mounted on angled stainless steel rods of 3 
mm diameter and were lowered into the tank between the edge 
of the model plate and side of the tank. Using clamps and a 
lead weight anchor located next to the tank, the hydrophones 
could be positioned anywhere in the tank. | 

Mounted to the source with the same size steel rod was 
an additional LC-10 acting as a receiver. With these two 
hydrophones colocated, the position of the source at the desired 
range and angle was accurately determined by measuring the 
travel time of the vertically propagating surface reflected signal 
and the travel time for the case of the rough surface or wedge. 
The receiver was positioned at the same depth as the source by 
measurement with a meter stick and the range was adjusted 
by measurement of the direct path travel time between source 
and receiver. A second hydrophone attached to the receiver 


would have allowed a more precise receiver depth 
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measurement, however, due to the probability of unwanted 
scattering, this method was not utilized. The positioning 
method utilized insured that the signal received was the 
specular scattered signal, 1.e., the angle of incidence equalled 
the angle of scatter. Once positioned, the source and receiver 
remained fixed throughout the data run. A set of reference 
data was then collected by reflection off the smooth water 
surface (reflection coefficient = -1). 

For the rough surface measurements, two additional data 
sets were required. First, to determine the plane surface 
acrylic plate reflection coefficient, the smooth plate section of 
the acrylic sheet was centered on the water surface between 
the source and receiver and reflected amplitudes in the specular 
(mirror) direction were recorded for a range of frequencies. 
second, to determine the rough surface reflection coefficient, 
the rubble surface of the model was centered by measuring the 
backscatter slant range as described above and the forward — 
specularly scattered pressure amplitudes were recorded. 

For measurements of the forward scattering from a single 
wedge fastened to the underside of a smooth acrylic plate, only 
one additional data set was required — that of the forward 
specularly scattered pressure off the wedge. The positioning of 
the wedge equidistant from source and receiver and at the 


desired grazing angle was accomplished as described above. 
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2. signal Processing 
For the effective frequency range e S0p 0 Kiera 

frequency resolution, Af, of at least 5 kHz was required in 
order to discern the peaks and troughs of the plotted data. The 
FFT algorithm used for processing the time domain pressure 
response required ‘power of 2' sample points. For the minimum 
frequency resolution of 5 kHz, the number of points used for 
the FFT was 256. Therefore, the actual frequency resolution, 
computed using Equation 5-1 was 3906.25 Hz. The 
transmitted signals were sinusoidal and were multiples of 


5906. Zon z 


sampling frequency 1,000,000 


= 3906.25 Hz (5-1) 
number of samples 256 


Due to a typical interference-free time window of less than 
100 usec, it was possible to collect only 64 time domain 
pressure sample points. Therefore, it was necessary to 
zero-pad the remaining 192 bins used for the FFT. 

The received signal was amplified 2000 times and filtered 
from 2 kHz to 300 kHz by the Ithaco Low Noise Preamplifier, 
then amplified another 20 dB and sent to the analog to digital 
converter and signal processor. Conversion to the frequency 
domain for analysis was performed on the signal average of 300 


time domaln measurements. 


60 


Pee oULITS AND DISCUSSION 
1. Results of Reflection from a Level Ice Model 
The plane wave reflection coefficient for a smooth surface 
is defined as R = p / Pin where p is the reflected pressure from 
the smooth acrylic plate and pj, 1s the incident pressure. The 
term p;, can be determined by the perfectly reflected pressure 
from the water-alr interface. 


For spherically diverging waves we use 


p (R, + R,) 


R, = (5-2) 
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In the above equation, 


R, = Slant range from source to specularly reflecting 
point on the plate 


R, = Slant range from plate to receiver 
p = Reflected pressure. The pressure of the received 
signal reflected trom the water-smooth plate 


interface as measured at range Ro. 


Py = Reference pressure. The pressure of the received 
signal reflected from the water - air interface. 


R, = Range at which reference pressure is measured. 
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As described in the procedure above, for each grazing angle, 
Rl = R2=0.5R,. Therefore, the range dependence drops out 
and the reflection coefficient becomes the ratio of the scattered 
received pressure divided by the reference reflected pressure. 
The reference pressure for both cases is the reflected pressure 
from a smooth water surface. | 

Measurements of the smooth plate reflection coefficient 
were made at five different grazing angles between 30 to 70 
degrees. For each angle, pressure measurements were collected 
at frequencies of 30.3, 50.8 and 74.2 kHz and are tabulated in 
Table II]. The reflection coefficients computed using Equation 
5-2 are plotted as dark squares for the above frequencies in 
Figures 5.1, 5.2 and 5.3, respectively. The solid line is each 
plot is the predicted pressure reflection coefficient from 
McCammon and McDaniel's computer model of reflection from 
arctic plates using as inputs the parameters of the bulk physical 
constants of the laboratory acrylic plate. Those parameters had 
previously been measured in our laboratory samples of acrylic 
material. The computer model predicts a single reflection null 
which is independent of frequency at 47° angle of incidence. 
The measured reflection minima range from 44° to 56° grazing 
angle. Since acrylic is a homogeneous viscoelastic material, the 
magnitude of the reflection coefficient at the null will decrease 


with increasing frequency. There is some evidence of this. 
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TABLE III 


SMOOTH PLATE REFLECTION COEFFICIENT 


REPLECTION CORPFICIENT 


Frequency (kHz) e005 50.8 74.2 
Grazing Angle 
(degrees 

oo OS io) 10 
a4 Ono D795 Gee> 
49 One Ce 0.80 
56 0 0.68 180 
63 heh eo) 0.94 
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On the whole, the laboratory results using a point 
source compare well with the computer model predictions for 
plane waves. Better agreement between theory and experiment 
could have been achieved by adjustments of our input speeds 
and attenuations. For example, a smaller compressional 
attenuation would yield deeper troughs in the computer 
calculations. We preferred not to make these adjustments. 
Although it was not our Intent to verify the accuracy of 
McCammon and McDaniel's computer model, our results do 


Show good agreement with their model predictions. 


2. Results of Coherent Soecular Scatter from a Model of a 
Rubble Field 


The plane wave coherent reflection coefficient for a 
randomly rough surface is defined as Ry = <p>/p;, where <p> is 
the average pressure of the forward scattered received signal 
measured for many realizations of the statistically rough 
surface. This concept is a generalization of the smooth plane 
surface reflection coefficient, R., which is the reflection 
coefficient derived frorm p, the pressure of the received signal 
reflected from a water -— smooth acrylic plate interface. 


Equation 5-3 defines R; for spherically diverging waves. 


R- = <p> —————— Ce 


R4, Eo Po and Roy are as defined for Equation 5-2. 

Figure 5.4 1s a plot of the mean of the square of the ratio of 
the specularly—scattered pressure relative to smooth plate 
reflected pressure, i.e., <(R-/R.) 42>, for different frequencies 
and angles of incidence. The abscissa of the plot is the surface 
roughness parameter, g, for specular scatter. This parameter is 
a dimensionless quantity which incorporates the rms surface 
height, sound frequency and incidence angle, 9;. Instead of 
scattering from different surfaces in order to get the average 
reflection coefficient, averaging has been performed over a range 
of angles and frequencies for bands + 0.25 g. The data, before 
averaging, are listed in Table IV. 

The measured data is plotted in Figure 5.4 against two 
theoretical curves: the solid curve is a plot of Equation 4-10 for 
a Gaussian PDF of heights and assuming R, is equal to 1. We 
have used the measured rms height deviation of the roughness 
elements for the assumed gaussian standard deviation. The 
dashed curve is a plot of the numerical integration of Equation 
4-8, Incorporating the measured pdf of the rough surface and 
assuming R, = 1. There is little difference between these two 
predictions which suggests that the non-Gaussian PDF of the 
model (and of the Arctic) does not cause a large error in the 
predicted coherent specular scatter. Analysis of the plot 
indicates that the Helmholtz theory provides reasonably accurate 
predictions of the specularly scattered rms intensities for values 


of g <.1, 1.e., low roughness. However at these grazing angles 
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Fig. 5.4 Specular scatter at model of water-rubble ice 

bottom interface. Smooth curve represents the theoretical 
value for the coherent component of scatter from a Gaussian 
surface. Dashed curve represents the theoretical value for the 
coherent component of scatter from the measured statistical 
surface. In both theoretical cases, Re is assumed 1. Dots 


represent mean of experimental results. Bars indicate standard 
deviation measured over several frequencies and angles. 


Boldface number indicates number of data points averaged. 
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dB Loss per Bounce 


MEAS IES PV 
PRESSURE RATIOS 
Grazing raed Smooth Plate/ Rough Plate/ Roughness 


Angle (KHz) Air to Water Air 10 Waterm Parameter 
Interface, <Rs> Interface, Rr (g) 
39 degrees 
SiS 0.98 0.80 On 71S 
Ca 0.88 Oreo 0.269 
Soe 0.94 0.44 0,352 
43.0 1.00 Osa 0,402 
46.9 1.00 0.69 0.478 
\Uec Oo 0.62 07561, 
54.7 1.00 C257, 0.650 
58.6 1.00 Un 5 0.746 
O25 i 00 Oey 0.848 
66.4 1.00 0.72 0.958 
70.3 1,00 0,60 Os 
74.2 1200 Oot lel OG 
Jom 1.00 0.12 1-925 
82.0 1.00 0.06 1.460 
49 degrees 
' gies 0.98 0.86 0.306 
Son 0.97 0.86 0.387 
39.1 0.96 ysis 0.478 
43.0 0.95 0.75 C573 
46.9 0.93 0575 0.687 
50.8 0.92 0.70 0,806 
54.7 0.89 O72 07935 
5870 0.89 Oey 2 1307S 
6215 0.86 0.76 1 228 
66.4 0.76 0.87 Loe 
LORS 0,82 0.79 1.544 
74,2 0.80 oo 76 172 
Tou Cay Qaey La77 
82.0 Ce 0.48 2, 100 


70 


aa eee eee” 


TABLE IV (cont'd) 
PRESSURE RATIOS 


Grazing ieee Smooth Plate/ Rough Plate/ Roughness 
i 


<— La i — ee —— ee eee ee ee eee eee 


Angle (KH Air to Water Air to Water Parameter 
Interface, <Rs> Interface, Rr (g 

63 degrees 
Ses iCe 0.88 0.427 
6572 0.90 0.89 0.539 
39.1 1,00 Ors 0.666 
43.0 0.91 0.85 C2305 
46.9 0.93 0.79 0.958 
50.6 Ons 0.70 1.124 
54.7 0.89 0.62 0S 
Soro Oe O52 1.495 
62.5 0.88 0.44 7 OL 
66.4 0.95 0.40 1.920 
7055 0.87 0.40 Z 152 
74,2 0.94 0.44 SM 
ol 0.89 0.46 Bs So 
82.0 0.89 0.40 2A 


PA 





for g > 1 there 1s a noticeable increase in the measured 
reflection coefficient relative to both theories for the coherent 
component alone. This trend is due to the increasing 
contribution of incoherent scattering from the steep-sloped 
roughness elements. The relatively large variance in the 
measurments for 1 < g < 2 15 due to the sensitive dependence 


on angle near the resonance. 


3. Results of Forward Specular Diffraction from a Model of a 
Ridge Keel 


Forward scatter from a single ridge keel mounted on a 
viscoelasctic smooth plate requires consideration of both the 
modified Biot-Tolstoy theory for forward diffraction from the 
ridge keel and the McCammon and McDaniel computer model 
predictions for reflection from a smooth viscoelastic plate. For a 
plate with a ridge mounted to its underside, one possible method 
to compute the energy reflection coefficient 1s defined by 
Equation 5-4. This method weighs the contribution of each 
theory's results by the percent of the total ensonified area 


which applies to the theory. 


Energy = (P+ wade) (A...) a (Peas out) (Ap) 
Reflection rn (5.4) 
Coefficient G25) (A,) 


In Equation 5-4 it Is assummed that: 


1. Ay = total scattering area. 


/2 


2. Ay, = area of wedge ensonified. 
5. Ap = area of smooth plate ensonified = Ay -— Ay. 
4. There is incoherent energy addition since: 


a. Wedge scatter is independent of plate scatter. 


b. The region of scatter is oe at the 
49th Fresnel Zone, at which region the zone 
macht is alcOulls Ome chi. 


The size of the ensonified elliptical area was determined by 
computing the semi-major and semi-minor axes based on the 
geometry of the experimental set-up and length of the time 
window. 

In this experiment all three edges of the ridge keel are 
ensonified by the incident sound energy, therefore, each of the 
three edges of the ridge keel must be considered when 
computing the total forward diffraction. The forward 
diffraction measurements are reported in terms of loss per 
bounce, i.e., 20 log (received pressure amplitude of sound 
diffracting from the ridge keel / received pressure amplitude of 
sound reflecting off a smooth water surface). Table V contains 
the measured data used to generate the above three plots. 
Figures 5.5, 5.6 and 5.7 plot measured results and compare the - 
results with the theoretical predictions for grazing angles of 27°, 
45° and 57°, respectively. In the upper half of each figure 1s 
an insert showing a source emanating acoustic energy which 1s 
being diffracted by each of the three edges (discontinuities) of 


the finite ridge keel. Reflected diffractions have not been 
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TABLE V 
WEDGE FORWARD SCATTER DATA 


Frequency(kHz Loss per bounce (dB re perfect reflection) 
Grazing Angle(degrees) 
33 45 63 

Ss ere S = oe =a 
Sone = ico =— -4.50 
Soa 2165 -4.30 -5. 64 
43.0 -1.94 =—4 2 -4.31 
46.9 rou -6.45 =o Oo 
S0ec -2.24 —6, © =GRa9 
54.7 —2.08 —2 6 oe 
58.6 -1.64 5 510 eS 
62.0 =17 16 -4,88 =. 
66.4 Sno =o oul ous 
70.3 = -3.13 -3. 56 
74.2 0.00 = eae -—4.38 
Kor meee -2.94 —Oueo 
82.0 =O =5 05 =O 15S 

Mean loss at each ie oor 4.19 

grazing angle aU 6G eS Sle aalee 
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Figure 5.5, Temporal and spectral forward scatter from the laboratory 
Wedge. Solid line is theoretical rigid wedge diffraction from the wedge keel. 
Dashed line is the theoretical combination of B-T diffraction and M&M 
reflection. Circles are experimental results. Grazing angle is 33 degrees. 
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Figure 5.6. Temporal and spectral forward scatter from the laboratory 
wedge. Solid line is theoretical rigid wedge diffraction from the wedge keel. 
Dashed line is the theoretical combination of B-T diffraction and M&M 

reflection. Circles are experimental results. Grazing angle is 45 degrees. 
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Figure 5.7. Temporal and spectral forward scatter from the laboratory 
wedge. Solid line is theoretical rigid wedge diffraction from the wedge keel. 
Dashed line is the theoretical combination of B-T diffraction and M&M 

reflection. Circles are experimental results. Grazing angle is 63 degrees. 


Ty 


considered here. The upper plot is a Biot-Tolstoy temporal 
description of the specularly scattered signal from a rigid wedge. 
The solid line on the lower plot shows the results of 
transforming the time data of the upper plot into the frequency 
domain. The dashed line near the top of the lower plot is a 
second theoretical curve which combines the Biot-Tolstoy 
predicted diffracted energy from the ridge keel and the 
McCammon and McDaniel predicted reflected energy off the 
ensonified area of the smooth plate, as explained in Chapter 4. 
The circles are our measured forward scattering results. These 
figures illustrate clearly that for this case of forward specular 
scatter from a single, finite wedge, the rigid wedge diffraction 
theory does not account for the total received acoustic energy. 
Also, a solution to the problem is not as simple as assuming the 
contribution of each component, 1.e., the ridge keel and the 
smooth plate, depending on its percentage of ensonified area. 

Analysis of Table V reveals that low frequency losses are 
only slightly dependent on frequency (within + 1.7 dB). There 
is a steady decrease in attenuation as the grazing angle is 
decreased. The difference between diffraction theory and 
measured results narrows as the grazing angle increases, as 
expected, since the ridge keel fills a larger percentage of the 
ensonified area as the grazing angle increases. 

One may summarize Figures 5.5 thru 5.7 by observing 
that the empirical results suggest that the loss per bounce at 


50 to 140 Hz in the Arctic decreases from about 4 dB at larger 
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grazing angles to about 2 dB at a 33° grazing angle. Also, the 
Biot—Tolstoy rigid wedge theory combined with a weighted 


reflection coefficient does not explain the measured data. 
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VI. BACKSCATTER 


A. PROCEDURE 

ocattered acoustic energy that propagates in the opposite 
direction from the incident acoustic energy is termed backscatter. 
To measure backscatter requires both the source and receiver to 
hie along the line connecting the source, receiver, and scattering 
obstacle. An LC~-10 transducer was used as the source to 
generate acoustic signals consisting of pulsed sinusoids. The 
propagating acoustic enerzy was scattered from one of the acrylic 
models described in Chapter III with the scattered acoustic energy 
measured by another LC-10 transducer located adjacent to the 
source transducer. Figure 6.1 shows the geometry used during 
the backscattering measurements. Typical measurements were 
made with the distance between the combination source/receiver | 
and the scattering obstacle of 50 cm. However, at steep grazing 


angles (angles greater than 50°) interference from the surface 


reflected acoustic energy required increasing this distance to 1 m. 
At the start of a new data run it was at times difficult to 
discern the backscattered signal from other scattering and 
reflections occurring in the tank. When this was the case, an 
aluminum sphere was positioned adjacent to the edge of interest 
and the resultant enlarged signal return was then used to isolate 
the signal of interest. Figure 6.2 is the time domain response 


from backscattering off the edge of an acrylic sheet. 
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Figure 6.1: Location of source and receiver with respect 
to the acrylic plate and scattering surface for the 
Beekscatter measurernents. 
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In the case of the open lead and ridge keel it was. impossible 
to isolate the response from individual edges with the available 
equipment. In those cases the source pulse length was 
maintained at approximately 150 ups. This insured an 80 us time 
window for the receiver in which the received signal was the 
sum of the diffracted signals fromm each of the individual edges. 
From this 80 ps window, 64 points were sampled at a frequency 
of 1 MHz and zero padded with 192 points to give sufficient 
frequency resolution after performing the discrete Fourier 
transform (DFT). The source frequencies selected in these 
measurements coincided with the bin frequencies from the DFT to 
minimize Gibb's effect (Clay and Medwin, 1977). The diffracted 
pressure from the sequence of edges was proportional to the 
magnitude of the response from the DFT. This was compared 
with the direct free field response over a fixed distance. For 
consistency, the same parameters were used for the backscatter 


measurements from the model of the rubble field. 


B. RESULTS OF BACKSCATIERUMEAcURG Vat iis 
All backscattering results are reported in terms of the 
backscattering strength (BSS) as defined in Medwin and Novarini 


(1981, p. 141). 
BSS = 10 log, [Geen eRe (6-1) 


where N is the degree of geometrical divergence, 0 < N 
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Figure 6.2: Time domain signal response for back— 
scattering off the lagging edge of an acrylic plate. Source 
and receiver are located 50 cm from the edge ata 
prazing angle of 30°. 
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is the acoustic intensity due to backscatter at the receiver; R is 
the distance between the scattering surface and the receiver, |, 
is the reference intensity at distance R,; and A is the area from 
which the scattering 1s measured. 

1. Spherical versus Cylindrical Divergence 

Prior to calculating the BSS it was necessary to determine 
if the divergence of the acoustic energy from a particular feature 
occurred either spherically or cylindrically at the range of the 
BSS measurement. The result affects the value of the term, N, 
in Equation 6-1 and hence the value of the backscattering 
strength. 

If the scattering occurs from a linear feature such as a 
lead or the keel of a pressure ridge, the divergence will be 
cylindrical if the pressure is measured close to the line scatterer 
and N will take on the value of 1. On the other hand, a 
randomly rough surface, containing individual point scatterers 
will result in spherical divergence, and N would take on the 
value of 2. 

Measurements were taken using the ridge model as the 
scattering feature with the source and receiver positioned as 
shown in Figure 6.3. The position of the source and the 
scattering angle to the receiver were kept constant while the 
recelver was moved to vary the range from the scattering 
surface. Hence, the pressure response of the receiver depends 
only on the range from the scattering surface to the receiver, 


Lee 
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Ridge Sail 





Paeure 6.5: Geornetry used in determining geometric 
divergence fromm a linear scattering element. R and ¢ 
were held constant at 0.5 m and 40° respectively, while 
rwes Varied. Dimensions of the ridge keel are as given 


in Figure 3.1. 
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on 
P< RR 


where R is the range from source to ridge, r is the range from 

ridge to recelver and n represents the type of divergence, i.e., n 
is 1 if the divergence is spherical and n is ’% if the divergence is 
cylindrical. Figure 6.4 shows a plot of log P vs log r with lines 

of slope ’ and 1 superimposed on the measured data. 

This figure clearly shows that, at the ranges used in the 
backscatter measurements, the divergence due to scattering from 
linear features is spherical, vice cylindrical as would be expected 
at shorter ranges. As described in Section [V.A, the impulse 
diffraction from a wedge is maximum at the point at which the 
least time path intersects the wedge; the contribution to the 
diffracted energy from points along the wedge on either side of 
the least time point decreases exponentially with distance from 
the least time point. Therefore, only a finite portion of the 
wedge effectively contributes to the diffracted energy which 
reaches the receiver. 

Using spherical divergence and converting intensities to 


pressures, Equation 6-1 becomes: 
BSS = 20 log, (eee k- ee ey) (6-2) 


For the measured data, P, was determined from the direct free 


field response of the receiver located at Rs = 1, or 2m. The 
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Figure 6.4: Results of the measurements using the 
geometry shown in Figure 6.3, where the range is r. 
Data indicates that the geometrical divergence is 
cylindrical («1//r ) at short range and spherical (««1/r) 
at longer ranges. 
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term A was taken as the entire area ensonified during the 64 us 
sampling period for the DFT. This was calcuated from 
information available in the averaged time domain response, that 


1s 
YA = SYnT(2t, + T) 


where c is the sound speed in the water, T is the record length 
for the DFT and t, is the elapsed time between triggering the 
source and the start of sampling for the DFT. 

Result kscatter from Simple Edges 

Simple edges consist of a single edge of an acrylic plate and 
can be divided into leading and lagging edges. A leading edge is 
defined when the source and receiver are located away from the 
plate, and 90° < @ < 180° as measured from the horizontal 
ice-water interface. A lagging edge is defined when the source 
and receiver are located under the plate, and 0 < 8 < 90°, 

Figures 6.5, 6.6, and 6.7 show the experimentally 
observed BSS plotted along with the theoretical BSS due to simple 
rigid wedge Biot-Tolstoy diffraction from the submerged leading 
edge. The time domain plot in these figures is the superposition 
of Equation 4-1 applied to diffraction from the submerged edge, 
diffraction from the submerged edge followed by reflection off the 
water-air interface, reflection off the water-air interface followed 
by diffraction from the edge, and reflection followed by diffraction 


followed by reflection. This gives three separate impulse 
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responses instead of four since Equation 4-1 gives identical results 
for diffraction-reflection and reflection-diffraction. The frequency 
domain plot compares the information in the computed time 
domain plot transformed to the frequency domain (by DFT) with 
the measured BSS. The sarnpling period for the DFT was chosen 
to give an integral number of bins (where one sampling period is 
a bin) between the least time arrivals from multiple or image 
wedges. 

For frequencies less than 50 kHz (83 Hz in the arctic) the 
overestimation of the BSS by Biot—-Tolstoy diffraction 1s due to the 
finite size of the plate edge, which at 8 mm is less than one 
quarter of the wavelength of the incident wave. Any 
discrepencies between the measured and theoretical data at 
frequencies greater than 50 kHz can be accounted for by mode 
conversion in the acrylic plate. Specular reflection does not 
contribute to backscatter. That is, some of the energy that 
would have been diffracted if the plate were rigid 1s converted 
into shear or flexural waves in the plate or refracted into the 
plate as a compressional wave; that energy does not appear in 
backscatter. 

In the case of the lagging edge, there will not be an image 
source. However, there will be some double diffraction, that 1s, 
diffraction around the edge with reflection from the water-air 
surface, followed by diffraction around the same edge in the 
direction of the receiver. This effect 1s expected to be small 


(Medwin et al, 1982) and will not be considered in this analysis. 
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Figure 6.8 shows the experimental backscatter from a 
lagging edge along with the predicted backscatter due to 
diffraction from the single visible edge. Only in the case of 15° 
grazing angle is there remote agreement between backscatter due 
to diffraction and the experimentally measured backscatter. At 
the other angles there Is significant backscatter from the edge 
that is not accounted for by rigid wedge diffraction alone. 
Possible mechanisms ce this additional backscatter could be: a) 
mode conversion of the incident acoustic energy into shear and 
flexural modes in the acrylic plate which are re-radiated into the 
water and b) compressional waves, internally refracting and 
reflecting in the acrylic plate and then re-radiating into the 
water. The discrepency is largest at grazing angles of 40 to 60° 
and is virtually non-existent at 15°. This behavior supports 
hypothesis b) above because the compressional wave in the 
acrylic plate approaches coincidence at the larger angles of 
Incidence, while the shear and flexural modes in acrylic do not 
achieve coincidence with the waterborne compressional wave. 
Even though coincidence is not achieved for the flexural and 
shear modes, some energy will be lost to these modes in the plate 
due to evanescent mode formation. This would account for the 
observed values of BSS that are less than those calculated for 
rigid boundary diffraction at the 15° grazing angle. 

3. Results from an Open Lead 
Backscatter measurements were taken from a simulated 


open lead with three separation distances and two grazing angles. 
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Figure 6.8: Temporal and spectral backscatter from a 
lagging edge of the acrylic plate. Lines are theoretical 
backscatter from a single edge at grazing angles of 1a, 30) 
45, and 60°. Points are the experimental data. 
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The lead was simulated by placing two acrylic sheets in the test 
tank and separating them by a fixed distance with the use of two 
balsa spacers. The separation distances used were 1, 2, and 

3 cm which correspond to Arctic dimensions of 6, 12, and 18 m, 
respectively. The two angles selected were a shallow grazing 
angle of approximately 15° and a steeper grazing angle of 
approximately 40°. The results of the measurements are shown 
Wemrisures 6.9 thru 6.14. 

Two fundamental concepts were used in evaluating this 
information: a) comparison with diffraction theory for diffraction 
from any visible edges and the water surface reflection of the 
diffracted energy, and b) comparison with diffraction theory plus 
the effect of multiple reflections in the water space between the 
plates. The effect of the multiple reflections was evaluated as 
diffraction of the incident impulse energy into the space between 
the plates where it is multiply reflected between the plates. To 
estimate the magnitude of this contribution, 50 percent of the 
remaining energy between the plates was assumed to be 
re-radiated back into the water to the receiver at discrete time 
intervals corresponding to the width of the plate separation. 

These were incorporated into the theoretical time domain response 
for diffraction as finite impulses whose magnitude included the 
effect of the initial diffraction and the effect of 50 percent leakage 
from the lead. 

The effect of diffraction from rigid edges alone is shown in 


the figures as the solid curves. The dashed curves include the 
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Figure 6.9: Temporal and spectral backscatter from an 
open lead with width 2 cm (12 m in the Arctic). Solid 
line fromm diffraction theory for two visible edges andeone 
surface reflection. Dashed line is diffraction theory plus 
the effect of multiple reflections in the lead. Points are 
measured values. Gragime angle sma. 
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Figure 6.11: Spectral backscatter fromm an open lead with 
width 3 cm (18 m in the Arctic). Solid line from 
diffraction theory for two visible edges and the surface 
reflection. Dashed line is diffraction theory plus the 
effect of multiple retlections tam ine me ao alee 
measured values. Grazing angle isu. = 
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Figure 6.13: Spectral backscatter from an open lead with 
width 1 cm (6 min the Arctic). Solid line is from 
diffraction theory for two visible edges and surface 
reflection. Dashed line 1s Giliracticone con meres eae 
effect of multiple reflections in the lead. Points are 
measured values. Grazing angle is 41°. 
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Figure 6.14: Spectral backscatter fromm an open lead with 
width 3 cm (18 m in the Arctic). Solid line is from 
diffraction theory for two visible edges and surface 
reflection. Dashed line is diffraction theory plus the effect 
of multiple reflections in the lead. Points are measured 
values. Grazing angle is 41°. 
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effect of multiple reflections in addition to diffraction. In 
computing the dashed curves a total of six reflections were 
considered. These results are consistent with those observed for 
the simple edges. Extending the results from simple edges at 15° 
grazing angle it 1s expected that BSS from both the leading and 
lagging edges that make up the lead will be predominantly from 
diffraction of the exposed edges. Due to the shallow grazing 
angle, coincidence with plate compressional modes which could 
result in enhanced BSS (from refraction and reflection) is not 
expected. However, mode conversion to evanescent 
compressional, flexural, and shear modes will occur resulting in 
reduced values of BSS. In addition, at the lower frequencies 
(<50 kHz) the observed BSS will be attenuated due to the finite 
size of the plate edge as discussed in the section on simple edges. 
At 41° grazing angle, conversion of waterborne compressional 
waves to plate compressional modes for the lagging edge of the 
lead with resulting refraction and reflection dominates the 
observed values of BSS. 
4, Results from a Ridge Keel 

Measurements of backscattering strength from a ridge keel 
were conducted in the same manner as those for the open lead. 
The same frequencies and approximate grazing angles were used 
so that a comparison between backscatter from the open lead and 
the ridge keel could be made. The results of the laboratory 
measurements and the analysis of the theoretical diffraction from 


a rigid wedge are given in Figures 6.15 through 6.18. 
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Figure 6.15: Temporal and spectral backscatter from the 
model of the’ ridge keel. The circles are measured back-— 
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Biot—Tolstoy diffraction theory. The grazing angle is 10°. 


Os 


Pressure Amplitude (arbitrary units) 


Boo da) 


Figure 
model 


800 


600 


400 


200 


~200 


-400 





~600 
-60 0 60 120 180 240 300 


Time (usec) 





-47 
30 35 40 45 50 55 60 65 70 75 80 85 


Laboratory Frequency (kHz) 


tH Ht 


50 60 70 80 90 100 110 120 130 140 
Arctic Frequency (Hz) 


6.16: Temporal and spectral backscatter from the 
of a ridge keel, Umerscoua neste enamel 


backscattering strength, lines are themed emia ecmien =m 
Biot-Tolstoy diffraction theory. The grazing angle is 29°. 


104 


1200 


800 


400 


-400 


-800 





Pressure Amplitude (arbitrary units) 
O 


=60 0 60 120 180 240 300 
Time (jtsec) 


BSS dp 





S035 401545 9S0 55° 60 65 70 75 80 85 
Laboratory Frequency (KHz) 


a ac 


SOC OMETCMMESONEESONNNIGON 110) (20 1350 140 
Arctic Frequency (Hz) 


Piguire ©)1/:)) lermporal and spectral backscatter frornm the 
model of a ridge keel. The diamonds are measured back-— 
scattering strength, lines are the calculated BSS using 

Biot-Tolstoy diffraction theory. The grazing angle is 47°. 


EUS 


seo 


250 


qZOo 





Pressure Amplitude (arbitrary units) 
© 


-500 
=o0 0 60 120 180 240 300 


Time (psec) 


BSS (dB) 





30 35 40 45 50 5S 60 65 70 75 80 85 
Laboratory Frequency (kHz) 


th ee re 
50 60 <0 80 = 90 100 110 120 130 £140 


Arctic Frequency (Hz) 


Figure 6.18: leémporal and spectral isaetseat ter 11 omamad= 
model of a ridge keel. The symbols are measured back- 
scattering strength, lines ere the calemiavcemibe. sins 
Biot-Tolstoy diffraction theory. The grazing angle is 63°. 


106 


At the shallow grazing angle of 15°, there is significantly 
more (10 dB more) BSS than can be accounted for by diffraction 
from rigid wedges. Due to the shape of the ridge keel, incident 
acoustic energy at the front edge of the ridge will be at a grazing 
angle of 42° (15° + 27°), which is approaching coincidence for 
compressional waves in the acrylic with respect to compressional 
waves in the water. These compressional waves can then be 
reflected and diffracted in the acrylic and re-radiated into the 
water resulting in the large BSS values observed. This same 
explanation also can account for the difference between theory 
and observation at 29° grazing angle (i.e. 56° grazing angle along 
front edge of the ridge keel). 

The results at 47° grazing =the: appear inconsistent with 
the results so far discussed. That is, if the observed BSS is 
dominated by compressional mode coincidence, then this should 
be the case for all frequencies since compressional waves are 
non-dispersive. Another effect that has not been discussed 1s for 
resonances to appear in the ridge keel. The dimensions of the 
ridge keel model are such that, depending on boundary 
conditions, the first resonance can occur at frequencies around 10 
kHz. Therefore the observed BSS can be affected by the 
resonances in the ridge keel and this can account for the 
fluctuations in the observed BSS. The appearance of resonances 
in the ridge keel model requires that the compressional wave 
energy be converted to compressional modes in the acrylic, since 


re-radiated flexural and shear modes are evanescent. 
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At the 63° grazing angle, the face of the ridge keel is 
normal to the incident acoustic energy, therefore there will be 
reflected energy in addition to diffraction, refraction, and 
re-radiation. Comparing these results with those at 47°, gives 
the contribution for reflection to BSS as approximately 15 dB. 

Figure 6.19 compares the observed BSS for the open lead 
and the ridge keel at shallow (15°) and moderate (41 to 47°) 
grazing angles. At shallow grazing angles the ridge keel 
contributes significantly more (15 dB) to BSS than the open lead. 
As discussed above, this is due to coincidence of compressional 
waves along the leading face of the ridge keel. This coincidence 
does not exist at 15° grazing angle for the smooth plate next to 
the open lead. At the moderate grazing angle, both the open 
lead and the ridge keel display comparable BSS. In this case, 
coincidence of compressional waves is present in both features. 

5. Results from the Rubble Field 

BSS measurements were made with the rubble field model 
and compared with actual Arctic BSS measurements. The 
measured BSS from the model of the rubble field for different 
grazing angles and frequencies is given in Table VI. In 
calculating the BSS, the area term in Equation 6-2 was taken to 
be just the area of the randomly rough surface that was 
ensonified, that is, the ensonified area of the smooth plate outside 
the rough area was not included. This differed from the 
calculations made for the open lead and ridge keel in which the 


area term included the smooth plate. This is justified since 
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Figure 6.19: Comparison of the measured BSS from the 
acrylic models of the open lead and ridge keel at shallow 
and moderate grazing angles. 
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smooth ice is an integral part of ridges and leads and its effect 
wil be included in actual arctic measurements, whereas 
backscatter from a large rubble area will not include smooth ice. 

The BSS was averaged over the laboratory frequencies of 
31 to 82 kHz (52 to 137 Hz scaled to the Arctic) for each grazing 
angle and plotted in Figure 6.20 along with data from actual 
Arctic measurements. The shallow grazing angle data are from 
Milne (1964) for frequencies of 25 to 100 Hz. These data were 
taken In a region of broken one year pack ice and therefore can 
be considered as being entirely from rough ice. The other set of 
data from Brown and Milne (1967) was taken in an area that 
contained 10 to 15% pressure ridges separated by rough ice. The 
low BSS at shallow grazing angles under these conditions can be 
attributed to shadowing by the ridge keels. Overall, these data 
Show the acrylic randomly rough surface to be a good model of 
an Arctic rubble field for low frequency acoustic backscatter. 

6. Composite BSS vs, Arctic Measurements 

The composite BSS for the acrylic model was calculated by 
converting the observed BSS for each feature (open lead, ridge 
keel, and rubble ice field) to backscattered pressure, weighting 
these pressures in accordance with the spatial occurrence of the 


ice features and converting the composite pressures to BSS. 
: BSS, /20 
BSS, = 20 log >, &,10°°™ 
i= 
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Figure 6.20: Backscattering strength as a function of 
grazing angle for the acrylic model of a rubble field 
compared with actual Arctic backscatter from rough ice. 
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feature, and &; is the spatial weighting factor. The values of € - 
used were: 0.05 for open leads, 0.125 for ridge keels, and 0.125 
for rubble ice fields. In this analysis, the BSS due to the smooth 
acrylic plate is taken to be zero. 

The composite BSS, scaled to five Arctic frequencies, 
calculated for a grazing angle of .15° is plotted in Figure 6.21 
along with BSS measured in the Arctic (Mellen and DiNapoli, 
1986) near 10° grazing angles. The plot shows excellent 
agreement between backscatter measurements made with the 
acrylic model and actual measurements made in the Arctic. The 
solid line included in Figure 6.21 1s from a recent computer 
model for backscattering from the Arctic ice (Green and Stokes, 
1985) based on diffractive scattering from small scale roughness 
on the flanks of randomly oriented triangular ridge keels. This 
model does not take mode conversion into account, which from 
the present work appears to be a significant contributor to 


backscatter from the Arctic ice. 


C. SUMMARY 

Based on the results of the backscatter measurements and the 
previous discussion, the following significant results are reported: 
(1) Mode conversion (predominantly refraction of compressional 
waves) significantly affects backscatter and can add or detract 
from backscatter due to rigid boundary diffraction by wedges and 
edges. (2) Backscatter from ridge keels at low frequencies and 


shallow grazing angles is more significant than backscatter from 
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Figure 6.21: Composite BSS computed for the acrylic model of 
the Arctic ice scaled to arctic frequencies plotted along with 

Arctic measurements of BSS. Curve is from a computer model of. 
arctic backscatter which does not consider mode conversion. 


open leads. The difference is attributed to coincidence between 
compressional waves in the water and compressional waves in the 
acrylic model along the relatively steep sloped ridge keel. At 
moderate grazing angles when coincidence of compressional modes 
is likely in both models, the observed backscattering results are 
comparable. (3) Laboratory backscatter data obtained from 
acrylic models, when weighted to account for statistical 
occurrence of smooth, ridged, open leads, and rubble ice yields 
backscattering strengths in good agreement with Arctic 


meaeasurements. 
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VIL. SMOOTH PLATE FLEXURAL WAVES 


The transfer of acoustic energy from the waterborne 
compressional wave to the ice medium in the form of ice plate 
flexural waves is studied to determine the significance of the 
energy loss of the incident waterborne compressional wave and 
the reradiation of flexural wave energy into the water boundary 
layer below the ice plate. For these measurements the smooth 
acrylic plate level ice model was used. The equipment set-up 
and signal processing was the same as for the previous 
experiments with the exception of the source and receiver as 


described below. 


A. PROCEDURE FOR THE MEASUREMENT OF THE SMOOTH 
PLATE FLEXURAL WAVE 


1. Measurement of Flexural Wave Speed 


The flexural wave speed was measured both with the 
plate suspended in air and floating in the water tank. In air 
the plate was laid on top of three 10 cm by 10 cm by 1.5 m 
wooden beams and measurements were collected in the area 
between the beams to minimize mass loading effects. For both 
measurements the plate was driven by an LC-32 transducer at 
25 volts. The transducer was held lengthwise in contact with 
the plate by resting a small bag of lead beads on top of the 
transducer. The purpose of this was to improve the driver 


system efficiency and prevent the transducer from Just 
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vibrating without putting energy into the plate. The receiver 
was an accelerometer. The accelerometer was screwed into a 
Pmewcin by 2.5 cm by i cm bDiece of plexiglass. Tackiwax® was 
used to stick the block of plexiglass to the plate. The signal was 
pulsed at three cycles per pulse with 100 msec pulse separation. 
Arrival times at 5 kHz intervals for a frequency range of 20 to 
80 kHz were collected for source to receiver separations of 50, 
45, and 60 cm. All data were stored on magnetic disk. 
2. Measurement of the Boundary Wave 

The smooth acrylic plate was floating in water for these 
measurements. The source was an LC-10 hydrophone positioned 
in the water below the plate. The receiver was an LC-10 
positioned below the plate in the water at a horizontal range of 
0.5 m from the source. The source and receiver length axes 
were parallel to the plate. The receiver was mounted using a 3 
mm stainless steel rod to a device that was able to adjust the 
depth at 1 mm increments. The initial measurement was 
taken with the source and receiver in contact with the plate 
from below. Once the flexural wave was identified, then the 
source and/or receiver were lowered at depth increments until 
the signal-to-noise ratio was too low for accurate 
measurements. Data were stored on magnetic disk. 

3. Measurement of Flexural Wave Attenuation 

These measurements were made with the plate floating 

in the water tank. The accelerometer could no longer be used as 


a receiver since the plate-to-receiver efficiency depended 


aly, 


heavily on the bond, and the previous method of using tackiwax 
for sticking the accelerometer to the plate could not be 
duplicated from range to range with a measurable precision. 
Therefore, measurements were made using a Bruel and Kjaer 
Yo in condenser microphone. This limited the frequency range 
since microphone frequency response roll-off occurred above 20 
kHz. The microphone, suspended from a metal rod spanning 
the plate, was slid across the plate at a height of 2 mm. Three 
data runs were conducted, one each at 10, 15 and 20 kHz. 

For the 10 and 1S kz data runs tive Garam coimts een mea 
spacing, were collected between 1 - 1.5 m Interval spacing. 
Below 1 m the time of arrival of the airborne sound wave 
generated by the source (travelling at 350 m/s) interfered with 
the plate flexural wave (travelling at speeds less that 650 m/s) 
at frequencies below 15 kHz. At 20 kHz the attenuation was 
too great for measurements to begin at 1 m; but due to the 
faster flexural wave speed (750 m/s) at 20 kHz, there was 
sufficient signal separation to collect data between 0.5 - 1.0 m. 
Since these measurements were collected in air in a noisy 
environment it was necessary to collect the data at night when 
ambient noise conditions were much lower. An HP-3561A 
Dynamic Signal Analyzer was used to collect and process the 
data. RMS averaging of 1000 samples points was used to 
compute the received signal level. Attenuation measurements 


and analysis were conducted by LT. Michael Browne, USN. 
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Peer LE AURAL WAVE MEASUREMENT RESULTS 


The condition necessary to achieve efficient mode conversion 
between compressional waves in the water and flexural waves 
in the plate is to match the horizontal component of the 
wavenumber in the plate. For both our acrylic plate and for 
Arctic ice, plate flexural speed is less than water compressional 
speed at all frequencies, therefore, the condition for efficient 
mode conversion is never satisfied. We were, however, able to 
generate flexural waves in the acrylic plate. 

1. Results of Flexural Wave Speed Measurements 

Once the flexural wave was identified the measurement 
of its speed was relatively easy. Identification required making 
several data runs while varying only the frequency and 
comparing the received pressure amplitude time plots. Due to 
the dispersive nature of flexural waves the received signal due 
to the flexural wave would arrive earlier as frequency was 
increased since flexural speed increased with increased 
frequency. Figure 7.1 is a plot of the acrylic plate flexural 
wave group speeds measured with the plate in air and floating 
on water. As expected, the speeds in water were slower than 
those in air due to the increased mass loading of the water on 
the plate. The frequency dependence of the flexural wave speed 
illustrates that for our frequency range of interest we are 
operating in the ‘low frequency’ regime. The flexural wave 


speeds of the rough surface acrylic plate were not measured but 
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Figure 7.1. Flexural wave phase speed as a function of 
frequency for a 3 mm thick smooth acrylic plate. Data was 
collected for tne plate in air and flogiing on water. 
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speeds slower than those measured for a smooth plate would be 
expected due again to the effect of mass loading. 

Figure 7.2 illustrates clearly that the acrylic plate 
accurately models arctic ice for the study of flexural wave 
propagation. The two lines marked ‘Ice are the theoretical 
flexural wave phase speeds plotted by substituting the reported 
physical constants of arctic ice listed in Table 3.1 into Equation 
4.17, the low frequency approximation, and assuming a 1.5 m 
ice plate thickness. The solid line marked ‘Acrylic’ is a similar 
theoretical plot for an assumed 1.5 m thick acrylic plate. The 
dots are laboratory measurements of flexural wave phase speed 
for the 3 mm thick acrylic plate in which the frequencies have 
been scaled by the appropriate thickness/wavelength ratio. 

2. Results of Boundary Wave Measurements 

As discussed in Chapter IV, any modes re-radiated into 
the water from a eeeorne flexural wave in a plate will be 
evanescent, that 1s, the energy will decay exponentially with 
distance from the plate instead of propagating as a harmonic 
wave. Figure 7.3 is a plot of the measured boundary wave 
pressure relative to the pressure of the direct, spherically 
diverging waterborne compressional wave plotted as a function 
of the ratio of the depth of the receiver below the plate to the 
wavelength of the waterborne compressional wave. The upper 
dashed line on the plot is the predicted level for the condition of 
the source in contact with the plate. The lower solid line is for 


the source at a depth of 6 mm. The circles and dark blocks 
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Figure 7.2. Theoretical flexural phase speed for an infinite 
flat plate of 1.5 mthickness. Ice and acrylic curves are 
calculated using values listed in Table 3.1. Dark circles iaige 
measurements collected on a 3 mm thick acrylic plate and 
scaled to 1.5 m thickness. The tzzzz area Detween the ice 
curves represents the range of theoretical arctic ice fien aig 


phase speeds. Arctic ice to acrylic frequency and thickness 
ratio is 1:600. 
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Fig. 7.3. Measured boundary wave pressure relative to the 
pressure of the direct water wave plotted as a function of the 
ratio of the depth of the receiver below the plate to the 
wavelength of the waterborne compressional wave. The upper 
dashed line is the predicted level for the condition of the 
source in contact with the plate. The lower solid line is for 
the source at a depth of 6mm. Thecircles and dark blocks 
are experimental results. 
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are experimental results. For both cases, the experimental 
results are in excellent agreement with the Tolstoy (1986) 
theory. If these results were extended to Arctic conditions, this 
graph suggests that for a 50 Hz signal from a source 3.5 m 
below the ice, the strength of the boundary wave at a receiver 
6 m beneath the ice and at a horizontal range of 300 m would 
be down no greater than 35 dB relative to the direct path 
signal. The term ‘no greater’ is used because plate flexural 
wave attenuation in ice is expected to be less than that 
measured in acrylic, whicn is the topic of the next section. 

3. Results of Acrylic Plate Attenuation Measurements 

The attenuation of the flexural wave in the acrylic plate 

was measured as 1.5 dB/m times the frequency in kHz. A 
value for Arctic ice flexural wave attenuation was not found in 


the available lterature. 
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ville COMSEUSIONS 


The results of these experiments provide insight into what 
contribution specific physical Arctic ice features make toward 
the gross underice “reflection coefficient” at low frequencies. It 
has been demonstrated that acrylic “ice sheets’ can acoustically 
model these ice features in the low frequency regime. This 
therefore permits one to conduct detailed studies of scattering 
from the Arctic ice canopy under controlled laboratory 
conditions. This is a major advantage over field measurements 
due to the highly variable nature of the ice floes and the 


difficultly in measuring the profiles of the under side of 3ea ice. 


A summary of the results of these experiments are: 


1. For a smooth acrylic plate, the pressure reflection 
coefficient 1s approximately 0.95 except near the 
coincidence frequency, where it may drop to 0.7. 


2. For a smooth acrylic plate, an evanescent flexural 
boundary wave exists, agrees with Tolstoy theory (1986), 
and 1s significant. 


CN 


For a pressure ridge keel, low frequency forward scatter 
losses are nearly independent of frequency, and decrease 
as grazing angle decreases. 


4. For a rubble field, a prediction of coherent specular 
scatter can be approximated ee aussian 
probability density function (PDF) of ice bottom roughness 
rather that as the actual pdf of the rough bottom. 
Losses up to 7 dB per bounce can be expected for large 
roughness. 


5. Mode conversion significantly affects backscatter and can 


add or detract from backscatter due to diffraction 
by ridge keels and ice edges. 
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Backscatter from ridge keels at low frequencies and 
shallow grazing angles is more significant than backscatter 
from open leads. 


Backscatter from the acrylic rubble field is in excellent 
agreement with scaled Arctic measurements from rough ice 
at shallow and moderate grazing angles. 


Laboratory backscatter data obtained from the acrylic 
models, when weighted to account for statistical occurrence 
of smooth, ridged, and rubble ice, and the presence of 
open leads yields backscattering strengths in good 
agreement with overall Arctic measurernents of 
backscatter. 
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